Inhomogeneous post-inflationary ACDM cosmology as a moduli space expansion 
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We model the large scale late time universe as a ACDM cosmology driven by cosmological constant 
and perfect dust fluid. Our aim is to find new solutions in the matter and A epoch consistent with 
inflationary initial conditions, namely that to the far past in the matter era the cosmology tends to 
a flat FLRW solution. We identify the moduli degrees of freedom that parametrize the flat A-dust 
FLRW solution and then promote these moduli to slowly varying functions of the spatial coordinates 
and show how to solve the Einstein equations in a comoving gradient expansion, controlled by the 
cosmological constant length scale. Our initial conditions ensure that the approximation remains 
under control to the far past of the matter era, and to the far future of A domination. The solution is 
fully non-perturbative in the amplitude of the metric deformation, and we explicitly construct it to 
fourth order in derivatives. A general A-dust universe dominated by A in the future is characterized 
by a 3-metric and a stress tensor (with positive trace) defined on the future conformal boundary. 
The new cosmologies with inflationary initial conditions are characterized only by the boundary 
3-metric, the stress tensor being locally determined entirely in terms of that metric. 
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I. INTRODUCTION 

At our relatively late epoch the matter content of the 
universe is well approximated by a cosmological constant, 
A, and dark matter [l], 0|. We live at an interesting 
time, approximately at A-matter equality, before which 
stretches a very long matter era, and to our future an 
infinite period of A domination. In this paper, as in our 
previous work Q, we take the view that it is interesting 
to focus on solutions to the Einstein equations which de- 
scribe the late time dynamics of our universe. We wish 
to use the natural late time scale, set by the cosmological 
constant, in order to separate the features of our cosmol- 
ogy that are relevant at late times from those that are 
irrelevant. We will model the A-matter era of a ACDM 
universe on large scales as a cosmology with cosmological 
constant and perfect dust fluid. Our aim is to find new 
tools to characterize and construct the most general in- 
homogeneous anistropic A-dust cosmologies relevant for 
describing our late time universe throughout the A and 
matter eras. 

Taking only A and dust matter gives an expanding cos- 
mology which inevitably has a dust dominated big bang 
to the past. For a A-dust solution to approximate our 
universe we must impose that the behaviour in the early 
dust era is consistent with inflationary initial conditions, 
namely that going back into the far past of the dust era 
the universe is increasingly homogeneous, isotropic and 
spatially flat. Said another way, we wish to construct 
A-dust cosmologies which tend to a flat FLRW solution 
in the far past of the dust era, and differ from it only by 
a growing deformation. We term such a solution a post- 
inflationary A-dust cosmology. We note that for a flat 
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A-dust FLRW cosmology the only length scale is given 
by the cosmological constant, and we will indeed see that 
the presence of A will be crucial to our construction of 
solutions. 

Cosmological perturbation theory ^a-|7| provides the 
conventional tool by which one constructs such deforma- 
tions of flat FLRW. Despite the success of perturbation 
theory we believe it is still interesting to pursue a differ- 
ent construction of inflationary ACDM cosmologies. In 
this paper, we detail a new construction which instead 
takes the form of a derivative expansion. The construc- 
tion can be regarded as a moduli space approximation, 
where the natural moduli of the flat FLRW A-dust cos- 
mology are identified, and then allowed to have slow spa- 
tial comoving variations controlled by the cosmological 
constant length scale. In cosmological perturbation the- 
ory, taking growing initial perturbations is crucial in or- 
der to consistently describe early times within the per- 
turbative approximation. Similarly in our approach, this 
initial condition is again important in allowing control 
to be retained in the early matter era within the gradi- 
ent expansion approximation. Whereas at late times in 
perturbation theory, growing perturbations may grow to 
the point of non-linearity, in our gradient expansion, the 
construction of the solution is fully non-perturbative in 
the metric deformation from FLRW and therefore this 
does not present a problem. 

Previous works have utilised the Hamilton- Jacobi 
framework to construct derivativel y e xpanded cosmolo- 
gies in the context of inflation [6j-[lC)j. The work pre- 
sented here differs from these approaches in two ways. 
Firstly, we are constructing solutions for the duration 
of A and matter dominance, rather than an inflationary 
epoch. Secondly and most importantly, we impose phys- 
ical inflationary initial conditions. 

In our previous work Q we have discussed how the 
general ACDM cosmology that is A dominated at late 
times can be characterized by its late time asymptotic 
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behaviour, in particular the future conformal boundary 
3-metric and stress tensor. We find that the cosmologies 
we construct here are indeed A dominated at late times, 
and are characterized precisely by their future conformal 
boundary 3-metric, with their stress tensor being fully, 
and perhaps surprisingly, locally determined in terms of 
this metric. 

The construction we use is analogous to the AdS-CFT 
discussion of emergent hydrodynamics from deformations 
of the homogeneous black brane jTl], ■ There, the bulk 
solutions are constructed as a derivative expansion con- 
trolled by the AdS length scale, and are shown to be 
characterized by a boundary stress tensor that describes 
a hydrodynamic theory. This derivative expansion may 
be thought of as a moduli space approximation, where 
the moduli give the physical degrees of freedom of the 
homogeneous black brane. The key insight in that set- 
ting is understanding the boundary conditions, namely 
that the asymptotic boundary geometry is unperturbed 
and requiring the future horizon remains regular. In our 
cosmological setting we may regard the boundary con- 
dition that at early times the deformation from FLRW 
contains no decaying component as being the analog of 
the horizon regularity condition. However, an important 
difference is that in our setting the boundary metric must 
be deformed in order to obtain interesting dynamics. 

The structure of this paper is as follows. In the next 
two sections we detail the Einstein equations and briefly 
review our previous work on characterizing general A 
dominated A-dust cosmologies from their late time be- 
haviour. We then discuss the nature of inflationary initial 
conditions, and review these explicitly using cosmological 
perturbation theory. In section [V] we detail the moduli 
space approximation for deformations of flat FLRW A- 
dust cosmologies, show how to impose inflationary initial 
conditions, and compute it explicitly up to fourth order in 
comoving spatial derivatives. In section[VI]we discuss the 
form of the boundary stress tensor that characterizes the 
solutions. Finally in section IVlII we briefly describe the 
use of Ricci flow to provide a natural smoothing flow in 
the space of post-inflationary A-dust cosmologies, which 
gives flat FLRW as a stable fixed point. 



II. THE EINSTEIN EQUATIONS 

We will take the matter content to be cosmological 
constant A and cold dark matter. We choose units so 
8ttGn = I- Since we will later be working in a deriva- 
tive expansion, controlled by the scale of the cosmological 
constant, we will be interested in large scales in the uni- 
verse where we expect the dark matter to be well modeled 
by pressureless perfect dust fluid. To begin we write the 
spacetime metric in a Gaussian normal form, 

ds 2 = -i + a 2 {y)g ij (x,y)dx i dxA (1) 



where x 1 for i = 1, 2, 3 give spatial coordinates on a con- 
stant y time surface and H 2 = A/3, and would give the 
Hubble rate for the de Sitter geometry corresponding to 
the cosmological constant A. Our time coordinate y is 
related to the usual proper time t as y ~ e~ Ht , so that 
y — > is the far future, and increasingly positive y cor- 
responds to the further back in the past. This choice 
of time coordinate simplifies discussion of the late time 
behaviour of our cosmologies, and conveniently compact- 
ifies future infinity to the conformal boundary y = pro- 
vided gtj is regular there. 

For our normal coordinates the vector d/dy is tan- 
gent to timclike geodesies. Hence we have the freedom 
to choose these constant y surfaces to comove with the 
dust fluid. Then the stress tensor takes the form, 

T w = J^yl ( A + P) > T vi = ' T ^ = ~ ^ fl2 3u 

and we do not have to introduce a velocity field for the 
dust, which is an important technical simplification in 
what follows. In particular it will allow us to straightfor- 
wardly identify the appropriate moduli of the flat FLRW 
solution. 

We may write the flat FLRW solution by taking 
gij(y, x) — Sij with the scale factor a(y) obeying, 

d 2d 1 / d \ 2 o let I 

h— — I =0, = - „ „ Pflrw 

a ya 2 \a J a ya bH z y z 

where ' = d/dy and p FLr (w denotes the spatially homo- 
geneous isotropic dust density. The solution is a 2 {y) = 

(l — !2_h-2 ?/ 3 ) 4 ^ 3 ; with p a constant of integration so that 
~ p y 3 +O(y 4 ) asymptotically at late times, y -> 0. 
However, by the trivial global scaling y — > Ay, x l — > Xx l , 
we may bring the scale factor to the form, 

a 2 (y) = (l-y 3 ) 4 / 3 (2) 

without loss of generality, and from now on we take a(y) 
to be defined in this way. 

We now perform a 1 + 3 decomposition of the Einstein 
equations. The Einstein equations then give a scalar and 
a vector equation with respect to the 3-metric gij , 

.. f 1 a\ . 1 . . w /i la\ I 

9 + + - )9 - ~9ij9 V +6 = -TT^P 

\ y a J 2 \a yaj H 2 y 2 

V'//,, V,// = 
together with the tensor equation, 

9ij + (™ + 3^ fa + ~gij (g mn g mn ~ g 2 ) 

-9i m 9j m + ~99ij + ~ = 0-(3) 

Indices are raised and lowered with respect to g^, Vj is 
the Levi-Civita connection of g^j, Rij is its Ricci tensor, 
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and we define g-y = d y gij, g = g^gij and likewise g\j = 
dy9ij> 9 = 9 l °9ij- The dust equation of motion is implied 
by these Einstein equations due to the contracted Bianchi 
identity. 

The scalar equation simply determines the dust fluid 
density. The vector equation is non-trivial. Let us denote 
the vector, 

$; = V 3 g i:j - 

so that the equation is <£>i = 0. The tensor equation Q 
implies, 

8 v $ i +(-- + 3?)* i = ~g$ i (4) 
V y a) 2 

and hence, as we would expect for a 1+3 decomposition, 
provided this constraint is satisfied on some constant y 
slice, it will always be satisfied. Hence we may regard it 
as constraining the space of solutions. 

III. REVIEW OF LATE TIME EXPANSION 

In this section we briefly review the late time charac- 
terization of A-dust cosmologies that become dominated 
by A [3]. For such solutions the asymptotic behaviour 
was deduced by Starobinsky in the context of inflation 
(l3| , and studied carefully more recently in [l4|, HH • At 
late times the general solution which is A dominated goes 
as, 

a(y) 2 9r 3 (y,x) = g tj + (Rij - ^SijRj + y 3 hij 

+y *aV+y 5 a?) + 0{y% (5) 

where (jij{x) is a freely specifiable non-degenerate Eu- 
clidean signature 3-metric which gives the geometry of 
the future conformal boundary y — 0. All tensors in 
the above expression are taken to live on the 3-geometry 
defined by <?.y , and indices are raised/lowered and deriva- 
tives taken with respect to this boundary metric, and R\j 
is the Ricci tensor of cjij . The tensor hij (x) is also general 
up to certain conditions. It is symmetric and an analysis 
of (U|) shows it must obey the constraint, V 3 'fty— Vji = 0. 
The dust density goes as, p = —3H 2 hy 3 + 0(y i ) 1 and 
hence we require that the trace h{x) < everywhere for 
physical dust densities. The tensors a„-J and a?) are 
given in 3]. 

The choice of 3-metric g~ij and symmetric tensor hij, 
with negative trace and obeying \7 J hij — Vj/i = then 
fully determine the solution at late times. All higher 
terms a[™' (x) are tensors determined covariantly in terms 
of this data (g, h). Thus we may take this pair (g, h) to 
characterize the A and matter eras of all ACDM cos- 
mologies that are asymptotically A dominated in the far 
future. The cosmologies we construct later are indeed in 
this class. 



From the tensor hij we may simply define Tij — hij 
cjijh, so that Tij is conserved in the sense of a stress tensor 
V J Tjj = 0. The positive dust density condition requires 
that its trace is positive, T > 0. The trace is zero in the 
case of vanishing dust density p = 0. 

In the context of dS/CFT and cosmological holography 
p^ - [23| . a naive analysis would suggest that Tij should 
play the role of the dual holographic stress tensor. Indeed 
setting p = so that the only matter is the cosmological 
constant, one can obtain a holographic renormalization 
gro up from analytic continuation of the AdS-CFT result 
[2J, [25| , and the traceless Ty would indeed give the dual 
stress tensor. Since dust leads to a non-vanishing trace 
of the stress tensor, if a holographic dual theory did exist 
for ACDM one would not expect it to be a conformal field 
theory. 



IV. INFLATIONARY INITIAL CONDITIONS 

The data (g, h) characterize the most general A-dust 
cosmology that is A dominated. However physically we 
are interested in A-dust cosmologies that have arisen from 
inflationary initial conditions, solutions that are deforma- 
tions of flat FLRW by only growing mode perturbations 
in the early matter era. Of course in our actual universe 
we believe the matter era is preceded by radiation and 
inflation, and we would expect some decaying perturba- 
tions in the early matter epoch. However, from the per- 
spective of a late time observer basing physics around the 
current cosmological constant scale, and characterizing 
the solution in terms of late time data (g, h), the matter 
era appears to be extremely long, unnaturally so from 
the perspective of the late time cosmological constant 
scale, and the early matter era appears to be very well 
approximated by a matter dominated big bang with only 
growing perturbations. The amount of decaying mode 
perturbation to be added in due to the fact that really 
there is no matter big bang, but a smooth transition to an 
earlier radiation epoch is vanishingly small when thought 
about in terms of the late time data, being suppressed 
parametrically by powers of the A-matter equality scale 
to the matter-radiation equality scale. It is this para- 
metric separation of scales that controls the decoupling 
of relevant dynamics in the A-matter era of our universe 
from that which is irrelevant, such as the details of the 
inflation and radiation era. 

Generic late time data (g, h) will correspond to a uni- 
verse which at early times is far from FLRW, in the sense 
that it will contain 'decaying' deformations that become 
stronger as one goes further back into the early mat- 
ter epoch. The task of this paper is therefore to under- 
stand how to characterize and construct the subset of 
cosmologies which are described by post-inflationary or 
growing deformations of FLRW such that the deforma- 
tions become smaller as one goes to earlier and earlier 
times in the matter era. We will give a construction 
that is non-perturbative in the metric deformation from 
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flat FLRW and instead takes the form of a derivative 
expansion. However, we have found it helpful to under- 
stand the conventional cosmological perturbation theory 
description of such post-inflationary A-dust cosmologies 
and this is reviewed and developed in the rest of this 
section. 

For small metric perturbations about flat FLRW, our 
metric choice ((T|) describes the synchronous gauge, where 
the residual gauge invariance has been fixed by requiring 
constant y slices to comove with the dust. Following 
the conventional approach [6, 7] we consider scalar, vec- 
tor and tensor perturbations with respect to the spatial 
metric. The perturbation is then written, 

9ij (x, y) = Sij + Sgfj (x, y) + 5gJ j (x, y) + Sgfj (x, y). 
Let us define the differential operator, 



A, 



2 a 
— +3- 

y a 



then the tensor perturbation, dgfj, 



traceless conditions d l Sglj = 5g 
and lowered with respect to the Euclidean metric , 
the evolution equation, 



obeys the transverse 
0, with indices raised 
and 



D y d y 5 9% = 



1 



a 2 (y)H 2 



AM- 



The scalar part is decomposed as Sgfj(x, y) — A(x, y)Sij + 
didjB(x,y), where the linearized Einstein equations im- 
ply, 



A(x,y) = (j>{x) 
D v d v B = 



l {y)H 2 



4>{x) 



for the scalar decaying and growing modes. 

We may analytically solve the scalar perturbation for 
all y which behaves at early times as the growing mode 
above; 



B(x) = c(x) 



(x)fy 2 (2+(l + y 3 ) 



.fi fig 3 ]) 



2H 2 



3(1 - y 



3U/3 



3jy 3 

l-y3 



(0) 



where 2-F1 is a hypergeometric function, and j is the con- 
stant, 



J 



2 8 / 3 y^r[5/6] 
3 5 / 2 r[i/s 



0.304.. 



(7) 



and c(x) is a constant of integration and we have picked 
a convenient normalisation for the mode. As for the vec- 
tor perturbations, the constant c(x) simply leads to a y 
independent infinitesimal diffeomorphism, and hence we 
will set c(x) = 0. 

However, for the tensor mode we know of no an- 
alytic solution for the combined A and matter eras 
and hence cannot analytically match this late time be- 
haviour to the early growing behaviour above. How- 
ever, performing a Fourier decomposition Sgfj(y,x) — 

J dk aij(ki)u(y,k 2 )e lHkmX , we may numerically solve a 
mode with wavenumber hi which obeys, 



DydyU = 



a 2 (y) 



(8) 



We now consider the late time behaviour of perturba- 
tions that grow in the early matter era. Firstly consider 
our growing scalar mode above. At late times this mode 
behaves as, 



where <j)(x) is an arbitrary function, but only of the spa- 
tial coordinates. The vector sector 5gYj is dynamically 
trivial, leading only to a y independent infinitesimal dif- 
feomorphism, and we will ignore it here. 

In the flat FLRW solution @ the big bang occurs at 
y = 1. Defining the new time, e = 1 — y, we may describe 
the early matter era by considering small e, where e be- 
haves as a proper time coordinate. Then a ~ e 2 / 3 , and 
the tensor and scalar modes behave as, 



Decaying 
Growing 



e 

B ~ const 



7*3 



S 9ij 



const. 



The scalar perturbation leads to a perturbation in the 
dust density, so p(y, x) = p^°\y) + Sp(y, x) where p^°\y) 
is the FLRW dust behaviour, and the perturbation 6p 
behaves as 



Decaying : 
Growing : 



5p_ 
p(°) 
5p_ 



1 



,2/3 



B 



1 



2H 2 



3 -? 3 7 

y 4>- 



2H 2 



0(y 4 ). 



The tensor modes are more complicated to treat and we 
can only numerically match the early time growing be- 
haviour since we lack an analytic solution. In terms of 
the Fourier decomposition, at late times one finds, 



0(y A ). 



We may determine the leading long wavelength behaviour 
for b(z) analytically by solving (|8]) order by order in z, 



c 2 z 2 + 0{z 3 ) 

whilst wc numerically determine c 2 as, 
c 2 ~ 0.09. 

We may then rewrite this in position space as, 



(9) 



SgZ(y,x) = 



d 2 

2H 2 



V 3 K-Jp)+0(jU*))<l>,.iU). 
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for some transverse traceless data $jj(x). 

We see that the scalar and tensor modes combine to 
describe a perturbation to the boundary metric, so that, 

9ij(x) = S i: j + ((j>8ij + . 

In fact, had we included the vector perturbation and the 
scalar field c(x) in ([6]), we would have also obtained the 
contribution dn£j) for some vector field £i(x), which is 
simply an infinitesimal difieomorphism. We observe from 
perturbation theory that the data hij is then completely 
determined by the perturbation to g as, 

hij(x) = -3 {Sij + {4>5ij + $y)) 

We now make two observations. Firstly we note that 
the general perturbation that grows in the early mat- 
ter era can be characterized by the perturbation of g, 
with the perturbation of h being completely determined 
by this. In a sense this is to be expected, since every 
perturbation wavenumber comes in both a growing and 
decaying time dependence at early times, we expect by 
choosing one behaviour we halve the data required to 
characterize the perturbation. 

Secondly, we observe that if we truncate h%j to two 
derivative order we obtain, to linear order in the pertur- 
bations, 

Shij(x) = -t&gu + ^ (W + + 0(d 4 ) 

= -\*9« - % (sRv - ~SR5^ + 0(9 4 ) 

where Scjij — <p 8ij + is the perturbation to the metric 
and SRij is the perturbation to its Ricci tensor i2y. 
Based on these observations, we might conjecture that 
the general solution which asymptotes in the far matter 
past to FLRW, as required for a post-inflationary uni- 
verse, is characterized by a general 3- metric g, with h 
being completely and covariantly determined by g as a 
derivative expansion with leading term, 

hij(x) = ~gij - \ Rij - ^Rg l3 ^j 

+0(V 2 R,R 2 ). (10) 

In the next section we show how to view the boundary 
metric as giving the moduli for the FLRW solution, and 
then allowing the metric to spatially vary gives a mod- 
uli space approximation for the general inhomogeneous 
anisotropic A-dust cosmology. We will then confirm that 
the prediction for the boundary above is correct. 

V. MODULI SPACE APPROXIMATION 

Let us now exhibit the moduli parameters of the fiat 
FLRW A-dust solution. Firstly we perform an arbitrary 



time independent global coordinate transform of the spa- 
tial coordinates, so that, 

bijdx l dx ) g~ijdx l dx J 

where g^ is a symmetric, non-degenerate tensor which 
is constant in time and space. Any such tensor can 
be obtained by a suitable choice of coordinate transform. 
The flat FLRW solution is then given by, gij(y,x) = gij. 
Our linear theory discussion above shows that the infla- 
tionary cosmologies should precisely be characterized by 
the boundary metric gij. Hence it suggests the constant 
components of the tensor g^ here are in fact the mod- 
uli of the flat FLRW solution, and that we can consider 
them to vary slowly spatially in order to develop a moduli 
space approximation. 

We know that a general late time A dominated solution 
for A-dust matter is characterized in terms of two tensors: 
a boundary metric g and a symmetric tensor h which 
obeys a conservation equation. Since our moduli appear 
to be precisely the boundary metric, we clearly must re- 
quire some boundary conditions in order to dynamically 
determine h. Following our discussion above, we take 
the boundary condition that the solution is asymptoti- 
cally flat FLRW to the far past, which can be thought of 
as the late time physical requirement that the universe 
had an inflationary origin. In practice we shall see shortly 
that this corresponds to a growing mode boundary con- 
dition on a linear differential operator, in analogy with 
the cosmological perturbation theory calculation. 

Let us characterize the spatial variation of g by taking 
the radius of curvature of the 3-metric g to be character- 
ized by the length L, so that, 

Let us define a perturbation parameter A = 1/(H L) 2 in 
which to perform our derivative expansion. We claim for 
varying moduli, the solution takes the form, 

9ij(y, x) = gtj(x) + g§\x, y) + gg\x, y) + 0(A 6 ) 

where g-™ '(x, y) ~ 0(X 2n ) is a tensor, living on the geom- 
etry defined by g^ and being a function of y, and being 
determined in terms of the moduli data g^ by a covariant 
expression in g^ involving (2n) spatial derivatives. We 
note that the derivative expansion is constrained to have 
only even powers by spatial parity symmetry. 

An important point is that the decomposition of the 
metric into a leading part, g^, that depends only on x 
and subleading parts depending on both x and y is not 
unique. We will find it convenient to define this splitting 
by requiring that g^ [x, 0) = for all 1 < n. This implies 
that the conformal boundary metric is simply given by 
cjij , and is not corrected by higher orders in the derivative 
expansion. 

We will now explicitly give the tensor equation and 
its solution for the trivial zeroth order, and the orders 
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determining g}j in the expansion. We then comment quadrature. Defining the source functions /CO(j/) 
on the structure of higher orders. Having determined 



and Slj\x) 
equation simply as, 

D y d v g$\y,x) = fW( V )S§\x). 

We note that these equations are simply tensor versions 
of the ones that occur for the scalar perturbations dis- 
cussed before. A particular solution, ffjj" can be ob- 
tained by quadrature, 

,d*).-. ^ _ c(D^y 2 (2 + (l + y 3 ) ^[1,1;!;^]) 



(1 2) 

the solutions g\-' , we will then compute the derivative 
expansion of the dust density, 

p(y, x) = pW (y) + pW (y, x) + p^ (y, x) + 0(A 6 ) 

where again p^"'(x,y) ~ 0(X 2n ) and depends on the so- 
lutions cl^ for m < n. In addition, we must also confirm 
that the constraint equation is satisfied order by order, 
and it is convenient to expand the constraint $j as, 



-gHj- (Rij — jgijR) we recast the first order 



(y,x) = A (sV^x) + *?\y,x) + 0(A 6 )) 



where the constraint (y,x) ~ 0(A 2 ™) and depends 
on the solution g^ 1 for m < n. 



A. Solution at trivial order: O(A ) 

Decomposing the equation (|3]) into orders in A one ob- 
tains at the trivial zeroth order, O(A ), simply the Fried- 
mann equation for a(y). The constraint equation is triv- 
ial at this order too. The dust density from the scalar 
equation is simply given by, 



P w (y) 



6H 2 y 2 



1 d 

ya 



2, ,3 



\2H z y 

(i-y 3 ) 2 ' 



We note that this is just the usual FLRW dust den- 
sity time dependence. One might think that this density 
should vary spatially in the moduli approximation, with 
large magnitude variations being allowed, which would 
seem to be at odds with the fact that the leading order 
behaviour is just that of the undeformed FLRW solution. 
However, we emphasize that we have chosen a co-moving 
coordinate system with the dust, and observe that on 
a comoving slice the dust density at leading order re- 
mains constant. This, of course, does not mean that the 
dust density is only perturbed in a small manner from 
FLRW since the constant y slice, with induced metric 
a 2 (y)cjij(x)/y 2 , can be strongly deformed relative to the 
FLRW one, which simply has a Euclidean geometry with 
induced metric a 2 (y)Sij/y . 



B. Solution at two derivative order: 0(A 2 ) 

The first non-trivial order is at two derivatives, 0(A 2 ), 
where one finds the tensor equation yields, 



Dyd y g^ = 



i 2 H 2 



Rii — 



1 



;9ij 



R 



where the right hand side is the Schouten tensor of the 
boundary metric. We note that due to the simple tensor 
structure, this is really a first order ordinary differen- 
tial equation for g^\y,x) and can by simply solved by 



g\;> i "(jj,x) = s% ) ( x y 



3(1 -y3)i/3 



where 2-F1 is a hypergeometric function. This particular 
solution at early times behaves as, 



2/3 



O(e) 



where as before e = 1 — y and j is the constant defined 
in ([7]), and at late times goes as, 



g§' p \y,x) = S$>(x) 1 ? + OW). 



(0/ 



The homogeneous solution is given by, 
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(y,x) = dijix) + 



3y 3 
1 - y 



^bij{x) 



with constants of integration a~ij (x) , bij (2) which we may 
regard as tensors living on the geometry defined by gtj, 
and behaves as, 



9%' h) (y,x) = bij(x)- c + faix) - 2b ij (x))+0(e) 



1 



at early times, and at late times goes as, 

g$ h \y,x) = Oij{x) + Zbytoy* + 0( y 6 ). 

We must now impose boundary conditions to solve the 
system. As stated, at early times we wish to impose infla- 
tionary initial conditions, namely that the perturbation 
has no decaying behaviour in the early matter epoch. 
Hence we choose the homogeneous data bij(x) to cancel 
the 1/e behaviour of the particular solution. The remain- 
ing data a~ij simply represents a change in the boundary 
metric gij(0,x) — §ij + a,ij + 0(A 4 ) at order A 2 . As dis- 
cussed above, we are free to choose the boundary metric 
to be given by gtj and not be perturbed at any order in 
the expansion. Hence we require djj = 0. Thus we have 
chosen, 

g§\y,x) = k™(y)S§ ) (x) 
with time dependence, 

' 3y 3 , y 2 (2 + (l + ^) 2 F 1 [l,|;f;^])' 



feCO(y) 



-J 



1-2/3 



3(1 -1/3)1/3 
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Since the equations here are of the form in the scalar 
perturbation theory above, it is no coincidence that we 
have found the same time dependence as in ([6]). 



j(2,o) 



The constraint equation at this order, $ 
trivial, and yields, 



(i) 



0, is non- 



a 2 H 2 



V 2 4 } " ViV^W - 2V ( ^) + -gij V m i>„ 



(1) _ 



(1) 



-2 i?. 



where indices are raised/ lowered with respect to g^. J ij ~ ilimy j y%j g yVmny \y I J U*3 

Since V^'S^ — V7.c(i) _ n +v,io ;mni;ao m no*romt 



0, this implies the constraint 



is indeed satisfied by our leading order solution. 

The dust density is also modified at leading order, and 
we can compute, 



from the scalar equation, and hence, 



In particular, we find that the perturbation to the FLRW 
zeroth order density by the first order goes at early times 
as, 



3 2 / 3 
40i? 2 



R(x)e 



2/3 



0(e 5/3 ) 



which has the same growing form as the in the lin- 
ear calculation. Had we not chosen the growing mode 
initial condition, and therefore had a perturbation in- 
volving some decaying behaviour, we would have found 
pW/p(°) ~ 1/e at early times, and hence the universe 
would not have approached FLRW in the far past. We 
see explicitly that for our choice of initial conditions the 
deformation g^\y, x) is a small perturbation to the lead- 
ing behaviour for all times, including during the early 
matter era, and hence this order of the gradient expan- 
sion gives a controlled approximation for all the A and 
matter epochs. 



C. Solution at four derivative order: 0(A 4 ) 

At four derivative order, 0(A 4 ), the equations are a 
little more complicated, and are given as, 



where indices are raised/lowered with respect to gij and 
v,i = V m <j£|] — V^ 1 ', which vanishes given the form of 
<7^ given above. This then yields the solution, 

= * <a - ) (v)s,? ,s) +* (?, ' 6) (»)^ ) 



where, 



?(2,o) 



1 

IF 



v 2 a 



r v 2 %. 



1 



— R gij + -RRij + —R mn Rrnngij — QRi" l Rmj 



£(2,6) = - 



[ ^ R2 9ij - g^-Rij - gR mn R mn gi 3 



H 4 
u m u 

it; It 



and the functions &( 2 > a / b ) solve, 
D ydy k^ = ^-, D y d v k^ = {m{y)f (11) 

such that they have no decaying behaviour at early times, 
so fc( 2 > a / b ) ~ O(e ) as e -> 0, and they do not lead to a 
deformation of g^, so that at late times £;( 2 ' a / fc ) ~ 0(y 3 ) 
for y — > 0. Note that unlike 0(A 2 ) the late time fall off 
goes as 0(y 3 ) rather than 0(y 2 ). This is simply due to 
the fact that we know the 0(y 2 ) late time behaviour from 
the late time expansion (|5j) only involves 2 derivatives. 
Hence order 0(A 4 ) and higher cannot contribute to it, 
provided the boundary metric is fixed to be gij. 

Whilst we cannot compute the functions fc( 2 ' a / b ) ana- 
lytically, numerically they are simple to determine. We 
give plots of the functions in figure [TJ They have the 
following behaviour at late times, 



fc (2 ' a/b) = c a/b y 3 + 0(y 4 ) 



(12) 



where c a and c b will shortly be given numerically. How- 
ever we can derive a relation between them. Combining 
the equations in (jlll) we find, 



D v 9 v9ij (V,x) = Sy*' a) (y,x) + S [ i-'"'{y lX ) 

where and are sources for the same linear differ- 
ential operator, D, as at lower order, and depend on the 
solution to the lower orders. Explicitly these are given 



(2,b), 



D y (fc( 2 '"> - 2k^ - 



0. 



(13) 



with a solution k^ 2 ^ - 2fc( 2 < Q ) - fcWfcW = 0, to which 
we may add the homogeneous mode of D y . Since we 
require that fc( 2,a ) and k^ 2 ' b ^ do not diverge at early times, 
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FIG. 1: The functions k {2 ' a \y) (solid) and k {2 ' b) {y) (dashed). 

and vanish at late times, we find that no addition of 
homogeneous mode to this particular solution is allowed. 
Expanding this solution at late times we find, 

c b = 2c a . (14) 

We find numerically 

c a ~ -0.184 (15) 

and confirm to the same precision that the relation (fT4")) 
holds. As a further check, we may study the lineari- 
sation of equation (|11[) . In a late time expansion, the 
coefficient of y 3 agrees with that obtained numerically in 
the tensor-perturbation calculation, ([9]), since c a — — 2c2 
to the precision obtained. 

(2) 

The constraint equation at this order, <&\ = 0, gives, 

-\&ig$ n )g {l)mn (16) 

with indices raised/lowered with g~ij, and must vanish for 
a solution. On our solution above, this translates into 
the condition, 

which as discussed previously is satisfied because it is the 
solution to equation (|13[) which is consistent with our 
chosen boundary conditions. The density at this order 
behaves as, 

„?>(„,*) = -*y (j» + (-i + 4) - 

-» (,w faM-HK)) 

and in fact we can explicitly compute the leading order 
early time behaviour analytically as, 

J2) o5/3 „• / 1 \ 

P 6 J IP. Smn 1 f=.2 \ 2/3 , ^/ 5/3n 



Again we see that for our choice of inflationary initial 
conditions this order of the gradient expansion gives a 
controlled approximation for the full duration of the A 
and matter eras. 



D. Solution at general order 

We believe the structure of the tensor equation at 
0(A 2 ™) order for n > 4 will be, 

D y d y g^(y,x)=J2p in ' a) (y)S (n ' a) (x) 

a 

where the sum of terms on the right hand side is a source 
for the linear problem with operator D. The spatial de- 
pendence of the source terms enters in the expressions 
g(n,a) which are composed from the solutions of the lower 
orders, S^ m ' b ^ for m < n, by tensor contractions, ab- 
stractions with Rij, and differentiation by V. Likewise 
the y dependence of the source is given by the functions 
p(n,a) wmc h are composed of expressions involving the 
f or m < n, the scale factor, and differentiation 

d/dy. 

The required solution of this tensor equation is then, 

a 

where the functions g(™' a ) solve, 

D y d y q {na) = p {n ' a) 

and the homogeneous modes are chosen so that q( n '°) ~ 
0(y 3 ) at late times (so the boundary metric gij is not 
perturbed), and there is no decaying behaviour ~ 1/e 
at early times as e — > 0, so we expect g("' a ) ~ O(e ), 
ensuring the gradient expansion gives a controlled ap- 
proximation at early times. 

A subtle point is whether such a solution obeys the 

(n) 

constraint $■ — 0. We have seen that at orders n = 1, 2 
this is satisfied, but it is not immediately obvious that 
we should have expected this, and we have seen that 
at order n = 2 the constraint is satisfied in a non-trivial 
way. We now provide an argument that for our boundary 
conditions the constraint should always hold. 

Let us assume that we have satisfied all the constraints 
up to order n — 1, so ^ n \y, x] = for all m < n, and 

(n) 

we have found the n-th order solution <?>„■ to the tensor 
equations, and that as claimed above, at early times the 
solutions behaves as g!-™' ~ O(e ) for all m < n since we 
have selected inflationary boundary conditions at each 
order. Now consider the equation On the right hand 
side g ~ 0(A 2 ), and hence if we consider order 0(A 2n+1 ) 
for this equation, we find the simple condition, 
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which implies <I > |™' ) must behave as, 
^ (n) (y,x) = h(x) 



(i-y 3 ) 2 



for some constants of integration bi(x). Now de- 
pends on both g!-^ and all previous orders g^ with 
to < n (for example, see equation ([T5|) ). The form of the 
constraint equation implies that provided g ™^ ~ O(e ) 

(n) 

for all rn < n at early times, then should remain 
finite as e — > 0, since there are no explicit divergences 
in the constraint equation, and so a divergence in 
could only arise from a divergence in g>j (y, x) for some 
to < n. However, we see from above that for bi{x) 7^ 

the constraint $^™' ) ~ 0(l/e 2 ) at early times. Hence we 
deduce that for our inflationary boundary conditions we 
must have bi(x) = and thus the constraint is satisfied. 
It is interesting that the constraint is not automatically 
satisfied, and it is an artifact of our inflationary choice 
of boundary condition that ensures this. We may regard 
the fact that the vector constraint equation gives no ob- 
struction to the moduli approximation at arbitrary order 
as confirmation that our choice of normal coordinates to 
the future boundary, together with the choice that these 
comove with the dust, is indeed consistent. 



VI. BOUNDARY STRESS TENSOR 

The linear theory led us to conjecture that cosmologies 
which in the past tend to FLRW, in the sense that the 
initial perturbation is a growing one, are characterized 
only by the metric g~ij, and the tensor hij is completely 
determined in terms of it. Using our moduli construc- 
tion we have explicitly constructed such cosmologies as a 
gradient expansion up to 4 derivative order. Calculating 
the boundary tensor hij to this order yields, 



-3^ 

Cg 



777 ( f '" 



1-72 u " ^ViV^R 



-0(A 6 ) 



t9ij 



-RRi 



23 
32 



where c a is defined in (1121) and given numerically in (1151) . 
From this, a corresponding expression for the boundary 



stress tensor may be obtained, Tj. 



jh. We note 



that this is indeed of the form predicted by the linear 
calculation in equation pH|) . We note that the 4 deriva- 
tive term above does indeed vanish on a scalar perturba- 
tion, as required for consistency with the linear analysis. 
From a holographic perspective we would say that the in- 
flationary initial condition corresponds to the dual QFT 
being in a particular state, associated to the metric g~ij 
on which it is defined. 



We may also re-express the tensor Ty by generating it 
through the variation of a Euclidean action, S[g 1 ^], 

- _ 1 SStfi] 
13 ~ V9 S3 1 * ' 

Generating Ty in this way guarantees that it is auto- 
matically transverse, provided <S , [g 4J ] is diffcomorphism 
invariant. We find the action which yields the deriva- 
tively expanded hij in (JTTJ) is given by, 



sis'- 



d 3 x^/§ 



H 4 



H 2 



R 



0(A 6 ) 



In the case of AdS-CFT, there has been considerable 
study of the moduli space deformations about the flat 
homogeneous AdS-Schwarzschild solution [n], [ill [26| . 
There, requiring the conformal boundary metric to be 
unperturbed and the horizon to be regular, one finds that 
the deformed solutions are characterized by the solutions 
of a dual fluid dynamics. We note that in this case, even 
for a trivial boundary metric there are interesting config- 
urations of the boundary stress tensor. By contrast, in 
our cosmological setting, the stress tensor and solutions 
become trivial for a trivial boundary metric cjij = 5ij. 
We believe this is a reflection of the fact that the dual 
holographic theory is Euclidean rather than Lorentzian, 
and hence would be expected to have a more restricted 
dynamics. 

We note that in the Hamilton- Jacobi formalism one 
also obtains a functional for a 3-metric , and indeed 

at fourth order in a derivative expansion, the same com- 
bination of curvature invariants is found 10]. It would 
be interesting to explore this connection further. 



VII. SMOOTHING POST-INFLATIONARY 
A-DUST COSMOLOGIES 

Since a post-inflationary A-dust cosmology is charac- 
terized by a Euclidean 3-metric, it is rather natural to 
consider flows in the space of solutions. A canonical ex- 
ample of such a flow is that given by the 3-dimensional 
Ricci flow, dgij/dr = —2Rij, where t is the auxiliary 
flow time (unrelated to real time). Since this acts only 
on the data g~ij which characterizes our post-inflationary 
A-dust cosmologies then by construction it gives a flow 
in the space of these cosmologies. Furthermore, it has 
three interesting properties. Firstly it provides a well 
posed parabolic flow. Secondly, we may interpret the 
flow from a starting 3-metric as a short distance smooth- 
ing procedure, where a flow for a time tq will result in 
a smoothing of a range of scales in the boundary metric 
(and hence cosmology) with characteristic scales shorter 
than ~ 1/y/fo- Thirdly the flat FLRW solution, g^ = 5^ 
is a (linearly stable) fixed point of the flow. 



10 



The use of Ricci flow for averaging cosm olog ical data is 
also considered by Buchert and Carfora |27H29| . However, 
an essential difference here is that for our solutions the 
data associated to the stress tensor is determined by the 
3-metric - a consequence of imposing our inflationary 
initial conditions. Thus we need only to smooth the 3- 
metric data and not the extrinsic curvature too. 



VIII. SUMMARY 

As an approximation to describing the A and mat- 
ter eras of our ACDM universe we have solved the A- 
dust Einstein equations using a moduli space expansion 
and impose initial conditions consistent with inflation, 
namely that in the early matter era, the cosmology tends 
to a flat FLRW solution. Our solutions are fully non- 
linear in the amplitude of the metric deformation away 
from FLRW. Furthermore, precisely because of our choice 
of initial condition they yield a controlled approximation 
from the earliest times in the matter epoch right through 
to the far A dominated future. 

An important simplifying step was the adoption of co- 
moving coordinates ([!}, for which the moduli are the 
components of the comoving FLRW spatial three-metric. 
This data is then promoted to vary on the comoving 
slices, and consequently the Einstein equations are no 
longer satisfied. Solutions may then obtained by intro- 
ducing time-dependent correction terms to the spatial 
three-metric and the dust density, order by order in spa- 
tial derivatives. Following this procedure we explicitly 
constructed solutions up to fourth order in comoving 
derivatives, imposing the inflationary initial conditions. 
We have also described the structure of the solution we 
expect at general order, and in particular we have demon- 
strated that the vector constraint equation gives no ob- 
struction. Our construction appears to reply on the pres- 
ence of the cosmological constant today which controls 
all scales in our solutions. More generally the unnatu- 
ral smallness of this late time scale allows a decoupling 



of the details of the early universe, such as the inflation 
and radiation era, from the dynamics of the late time 
universe, the matter and A eras. As remarked in the in- 
troduction, the moduli space approximation we use here 
is strongly reminiscent to the derivation of hydrodynam- 
ics from gravitational solutions in AdS/CFT. 

In a general late time characterisation of ACDM so- 
lutions, as discussed in [3], a A-dust universe dominated 
by A in the future is determined by data on its future 
conformal boundary, namely a 3-metric and stress ten- 
sor. We have found that our choice of inflationary initial 
conditions determines the stress tensor in terms of the 3- 
metric. Whilst such a reduction in data is to be expected 
when one moves from a general solution to one obeying 
certain initial conditions, it is perhaps a little surprising 
that the expression for the stress tensor is a local one in 
terms of the 3-metric. 

Whilst one might have thought that characterizing our 
cosmology by using 'future' data rather than initial data 
would be practically hopeless if one was interested in par- 
ticular initial conditions, such as inflationary ones, we 
have seen that at least on large scales this is not the case 
at all. It then remains an interesting question whether 
we may use the construction presented here to sharpen 
large scale observational constraints at late times using 
late time observations, for example using the luminosity- 
distance vs. redshift relation computed in [3J, even when 
we make the assumption of inflationary initial conditions. 
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